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ABSTRACT 

Coagulation of submicron-sized dust grains into porous aggregates is the initial step of dust evolution in 
protoplanetary disks. Recently, it has been pointed out that negative charging of dust in the weakly ionized 
disks could significantly slow down the coagulation process. In this paper, we apply the growth criteria obtained 
in Paper I to finding out a location ("frozen" zone) where the charging stalls dust growth at the fractal growth 
stage. For low-turbulence disks, we find that the frozen zone can cover the major part of the disks at a few to 
100 AU from the central star. The maximum mass of the aggregates is approximately 10~ 7 g at 1 AU and as 
small as a few monomer masses at 100 AU. Strong turbulence can significantly reduce the size of the frozen 
zone, but such turbulence will cause the fragmentation of macroscopic aggregates at later stages. We examine 
a possibility that complete freezeout of dust evolution in low-turbulence disks could be prevented by global 
transport of dust in the disks. Our simple estimation shows that global dust transport can lead to the supply of 
macroscopic aggregates and the removal of frozen aggregates on a timescale of 10 6 yr. This overturns the usual 
understanding that tiny dust particles get depleted on much shorter timescales unless collisional fragmentation 
is effective. The frozen zone together with global dust transport might explain "slow" (~ 10 6 yr) dust evolution 
suggested by infrared observation of T Tauri stars and by radioactive dating of chondrites. 
Subject headings: dust, extinction — planetary systems: formation — planetary systems: protoplanetary disks 



1. INTRODUCTION 

Growth of submicron-sized interstellar dust grains into 
kilometer-sized planetesimals is the first step towards the 
formation of terrestrial planets and th e cores of gas giant s 
in protoplanetary disks ilMizunol [19801 iPollack et al.l U996). 
The formation of planetesimals is believed to involve the fol- 
lowing stages. (1) Initially submicron-sized particles coagu- 
late into larger but highly porous, fractal aggregates through 
Brownian motion and differential settling towards the mid- 
plane of the disk dWurm & Bluml [19981 iBlum et all [19981; 
iKempfet al. fl999l) . (2) As the aggregates grow to "macro- 
scopic" (mm to cm) sizes, the collisional energy be comes 
high enough to cause the c ompaction of the aggregat es {Blum 
l2004tlSuyama et al.ll2008l: iPaszun & Dominikl2009l) . (3) The 
compaction cause the decrease in the gas drag force acting 
on the aggregates, allowing them to concentrate in the mid - 
plane of the disk dSafronovlll9"69T: IGoldreich & Ward [19731) . 
the center o f vortices (Barge & So mmeriall 19951) . and turbu- 
lent eddies ( Johanse n et al.l 120071) . (4) Planetesimals may 
form within such dense regions through gravita tional insta- 
bility (ISafronovlll969T: IGoldreich & WardN 19731) or through 
further collisional gr owth fWeidenschilling & Cuzzi 119931; 
IWeidenschillh^[T995l) . 

However, a number of obstacles, often called "barriers," 
have been pointed out against the above processes. As the col- 
lisional compaction proceeds, the motion of aggregates rela- 
tive to the gas becomes faster and faster due to radial drift 
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dWeidenschilring|ll977l) a nd random motion induced by tur- 
bulence dVolk et al.lll980l) . The colli sion velocity can excee d 
10 m s -1 even without turbulence dWeidenschilUnd 119771) . 
but it is uncertain whether such high-speed collisions lead 
to the sticking or fr agmentation of the aggregates (the "frag- 
mentation barrier"; iBlum & Wurmjl2008t IWada et all 120091 ; 
iTeiser & Wurml 120091: IGuttler et aLlboidl) . Furthermore, re- 
cent laboratory experiments suggest that the compaction it- 
self causes reduction o f sticking efficiency and bouncing 
of co lliding aggregates (Blum & Wurrril 120081: IGuttler et all 
2010). This may halt dust growth before the aggre gates reach 
the fr agmentation barrier (the "bouncing barrier"; Zso m et al.1 
2010). Even without fragmentation, the radial drift can re- 
sult in the lost of dust particles from the entire part of 
the disk unless their gro wth proceeds very rapidly (the "ra- 
dial drift barrier"; e.g.. Br auer et al.1 12008). Besides, tur- 
bulence does not necessarily promote the local concentra- 
tion of dust, because it can effi ciently diffuse small particles 
dWeidenschilling & Cuzzilll993l) . 

On the other hand, astronomical as well as meteoritic ob- 
servations suggest that some barriers against dust evolution 
may do exist in protoplanetary disks. Mid-infrared excess 
observed for classical T Tauri stars (e.g.. iFurlan et al.ll2006l) 
implies a certain amount of small dust grains retained in the 
inner parts of their circumstellar disks over a million years. 
This feature cannot be explained by sim ple coagulation theory 
assuming perfect sticking efficiency (Dullemond & Dominik 
120051) . In addition, radioisotope dating of the most primi- 
tive chondrites supports that the formation of chondrules be- 
gan at least a million years after the formation of the Solar 
nebula (e.g. iKita et al.l l2000. 2005). This also suggests that 
the dust growth p rocess towards planetesimals is "inefficient" 
dCuzzi et al. 2008). Thus, for better modeling of dust evolu- 
tion in protoplanetary disks, knowledge on what prolongs it is 
as important as knowledge on what promotes it. 

Recently, one of the authors has pointed out another kind 
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of b arriers, namely, an electrostatic barrier due to dust charg- 
ing (Okuzumi 2009, hereafter O09). Protoplanetary disks are 
expected to be weakly ionized by a various kinds of high- 
energ y sources, such as cosmic rays ([Umebavashi & Nakano 
1981) and X -rays from the central star (Glassg old et aljl997t 
Igea & Glassgold 1999). As is known in plasma physics 
(Shukla & Mamun 2002), dust particles negatively charge in 
an ionized gas because free electrons hit to the particles 
more frequently than ions. This "asymmetric" charging im- 
plies possible existence of an electrostatic barrier against dust 
growth, but this effect has been ignor ed in the context of dust 
growth in protoplanetary disks. IO09I estimated how strongly 
the charging affects the collisional cross section between ag- 
gregates properly taking into account the weak ionization of 
the disks. The result shows that the collisional cross section 
can be strongly suppressed before the collisional compaction 
becomes effective. This is in clear contrast to the previously 
known barriers which act after the o nset o f the compaction. 

However, the simple estimate by IO09I assumed that dust 
aggregates grow with a narrow size distribution. In reality, 
size distribution is determined as a result of the coagulation 
process, and it has been unclear how the distribution evolves 
when the charging is taken into account. The number of small 
aggregates is particularly important because it determines the 
ionization state of the gas and hence the charge state of all ag- 
gregates dO09l). To ad dress thi s issue, i n our preceding paper 
(Okuzumi et al. 201 1, hereafter Paper J), we have numerically 
simulated how the size distribution evolves in the presence of 
the electrostatic repulsion properly takin g into account ion- 
ization balance in the gas-dust mixture (O09) a nd porosity 
evolutio n due to low-velocity sticking (Ok uzumi et al .1120091 
hereafter OTS09). We find that the outcomes can be classified 
into three types: 

(a) Unimodal growth. If the electrostatic repulsion is neg- 
ligibly weak, aggregates simply grow with a relatively 
narrow size distribution (Figure Qj a)). We refer to this 
growth mode as the "unimodal growth." 

(b) Bimodal growth. If the electrostatic repulsion is strong 
but nonthermal motion (e.g., vertical sedimentation and 
turbulence) dominates aggregate collision, some ag- 
gregates stop growing at a certain size while the rest 
continue growing by colliding with each other (Fig- 
ure 02b))- We call this the "bimodal growth." In this 
mode, growing aggregates dominate the total dust mass, 
but their negative charges are suppressed by the non- 
growing ("frozen") aggregates. Interestingly, in some 
cases, the presence of the frozen aggregates is even re- 
quired for the larger aggregates to continue growing. 
The bimodal growth thus demonstrates the importance 
of the dust size distribution. 

(c) Total freezeout. If the electrostatic repulsion is strong 
and thermal (Brownian) motion dominates aggregate 
collision, all aggregates stop growing at a certain size 
(Figure [TJc)). We refer to this type of growth mode as 
the "total freezeout." 

We have also obtained a set of simple criteria for which of the 
outcomes is realized under given conditions. These criteria 
allow us to predict how the initial fractal growth proceeds at 
different locations in protoplanetary disks. 

In this paper, we use the growth criteria obtained in lPaper 3 
to map a region where local fractal dust growth ends up with 



(a) Unimodal Growth 



o 

o ° o 

o O 



(b) Bimodal Growth 



o ° 
o o 



°°o ° 



o 



o 

O O 
o 



o./o 

° o 
o (^^J o o 



(c) Total Freezeout 



FIG. 1. — Schematic illustration showing three outcomes of the collisional 
growth of fractal dust aggregates charging in a weakly ionized gas: (a) uni- 
modal growth, (b) bimodal growth, and (c) total freezeout. In the unimodal 
growth, dust aggregates grow with a relatively narrow size distribution. In the 
bimodal growth, a certain number of aggregates stop growing due to the elec- 
trostatic repulsion and the rest continue growing by colliding with each other. 
In the total freezeout, all aggregates stop growing with a nearly monodisperse 
distribution. See Paper I for details. 



the total freezeout, to which we will refer as the 'frozen zone." 
This is the first step towards comprehensive modeling of dust 
evolution in protoplanetary disks in cluding dust charg ing to- 
gether wit h collisional compaction dSuvama et alj|2008l) . ra- 
dial drift dWeidenschilling|ll977t iBrauer et alj|2008[). bounc- 
ing dZsom et al.l2010l) . and fragmentation dBrauer et alj2008fc 
iBirnstiel et al]l2009h. In Section 2, we briefly summarize the 
analysis done in Paper I and present the growth criteria. The 
protoplanetary disk model used in this paper is described in 
Section 3, and the main results are presented in Section 4. In 
Section 5, we discuss potential mechanisms that could pre- 
vent dust growth in the frozen zone from completely frozen. 
A summary of this paper is presented in Section 6. 

2. COAGULATION OF FRACTAL DUST AGGREGATES 
IN A WEAKLY IONIZED GAS 



In this section, we outline the analysis done in PaperJ and 
introduce several quantities to write down the growth criteria. 

We focus on the first stage of dust evolution in pro- 
toplanetary disks where submicron-sized dust particles 
("monomers") grow into fractal aggregates. For simplicity, 
we assume that the size of the monomers is equal and treat 
it as a free parameter. The assumption of fractal growth is 
valid only when the collision energy is low and both com- 
paction and fragmentation is negligible ( Domini k & Tielensl 
ll997tlSuyama et alJ l2008). The limitation of this assumption 
will be shown in Section 4. 1 . Note that most previous stud- 
ies on dust coagulation ignored fract al evolution and treated 
aggregates as compact spheres (e.g., iNakagawa et alj 119811 : 
Tanak a et alJl2005b IBrauer et alJl2008l) . However, as we will 
see in Section 5.1, it is critical to properly take into account 
fractal evolution when analyzing the electrostatic barrier. 

Furthermore, we assume that the dust growth proceeds lo- 
cally. In laminar disks, this assumption is valid as long as dust 
grows into fractal aggregates, because the timescales of verti- 
cal settling and radial drift are much longer than that of local 
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growth. In turbulent disks, dust can be globally transported 
on a short timescale by turbulent mixing, but we ignore this 
effect in Sections 2-4. Effects of turbulent mixing as well as 
radial drift will be discussed in Section 5.2. 

The collision velocity between aggregates is assumed to be 
driven by thermal (Brownian) motion and nonthermal "dif- 
ferential drift." Here, differential drift refers to dust motion 
whose relative velocity has of the form Aud = g\Tf t \ —Tfp\, 
where t/j and ry 2 are the stopping times of colliding aggre- 
gates and g is the effective acceleration inducing the relative 
motion. In protoplanetary disks, main sources of g for small 
aggregates are stellar gravity towards the midplane of the disk 
and turbulent motion of the ambient gas (see Section 3). 

The charging mechanism considered in this study is the 
captu ring of ionized gas particles in a weakly ionized gas 
(O09). An important parameter characterizing the gas ion- 
ization state is the ionization rate £, which is the probability 
that a molecule is ionized into an ion-electron pair per unit 
time. 

2.1. The Kinetic and Electrostatic Energies 

Collision of charged aggregates depends on the kinetic en- 
ergy of their relative motion and the electrostatic energy. The 
kinetic energy is the sum of the thermal energy (~ k B T, 
where k B is the Boltzmann constant and T is the temper- 
ature) and the energy associated with the differential drift, 
E D = [M\M 2 /{My +M 2 )](Au D ) 2 /2, where M\ and M 2 are the 
masses of colliding aggregates. The electrostatic energy is 
defined as Ee = Q\Q 2 /{a\ + a 2 ), where 012 and Q\.2 are the 
radii and charges of the aggregates. As shown in Paperll the 
collision probability is significantly suppressed if Ee is much 
larger than k B T and E D . Therefore, the outcome of dust evo- 
lution is determined by how these energies increase as aggre- 
gates grow. 

As shown in lPaper fl aggregates grow with a relatively nar- 
row size distribution until the electrostatic repulsion becomes 
significant. For this reason, it is useful to evaluate Ed and 
Ee assuming that all aggregates have the same mass at every 
moment. Following Paper I, we will call this the "monodis- 
perse approximation." Within the fractal growth regime, this 
is equivalent to treat the collision products as bal listic cluster- 
cluster aggregates (BCCA; e.g.. iMeakinl 119911) . A BCCA 
cluster is characterized by the monomer number N = M/mo 
and radius a sw ac,N [ / D , where mo and aq are the mass and 
radius of the monomers an d D ~ 1 .9 is the fractal d imen- 
sion of BCCA clusters (e.g.. lMukai eTafl ll9"92: OTS09]). We 
simply set D = 2 in this paper. Indeed, A^-body simulations 
suggest that aggregates have a fractal dimension close to 2 
unle ss highly unequal-siz ed collisions dominate their growth 
(see lOkuzumi et afl l2009). Thus, D ~ 2 is a good assumption 
for aggregates growing with a relatively narrow size distribu- 
tion (and without collisional compaction). 

The drift velocity Auo is given as follows. Focusing on 
early stages of dust evolution, we assume that the radii of dust 
aggregates are much smaller than the mean free path of gas 
molecules. Under this assumption, the stopping time ry of an 
aggregate can be given by Epstein's law, 
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where p g is the gas density, m g is the mass of gas molecules, 
and A is the projected area of the aggregate. Under the 
monodisperse approximation, the drift energy normalized by 



k B T can be written as (see Equation (37) of iPaper fl) 

p - Ed - t 2 ^ 3 ™ 
= ksT ~ ytfTV) 2 ' (2) 

where A = A j r Ka\ is the normalized mean projected area, e is 
the ratio of the standard deviation to the mean of the mass-to- 
area ratio for given M, and fo is the drift energy of monomers 
normalized by k B T . Denoting the bulk density of monomers 
b y pn = 3m n/47rfln, the definition of f D reads (Equation (28) 
of Paperll). 
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We deter mine A according to a fitting formula for BCCA 
clusters bv lMinato etafl (120061) . 



A(N) = 



12.5N a685 exp(-2.53/Af a092 °), N < 16, 
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N> 16. 



(4) 



Note that the mass-to-area ratio N/A approaches to a constant 
in the limit of N ^> 1. T his is a general feature of fractal ag- 
grega tes withD < 2 (see lMeakin & Donnll988l:lMeakin et alj 
1989). In Paper I, we have calculated e for numerically cre- 
ated BCCA clusters and found that e ~ 0.1 depending on N 
only weakly. We will assume e = 0.1 independently of N. 

The dust charge Q is calculated in the following way. We 
introduce a dimensionless, negative surface potential ^> = 
-Qe/aksT. I f the gas surrounding the dust is fully ionized, 
* is given by (Spitz a1ll94l1IShukla & Mamunll2002l) 



1 



1 + * 



■exp* = 0, 



(5) 



where m^e) is the mass of ions (electrons), and is their 
sticking probability onto dust surfaces. We write the solu- 
tion to Equation (|5]) as ^oo. Assuming m, = 24mn (corre- 
sponding to the mass of Mg + ), Sj = 1, and s e = 0.3, we obtain 
$oo~2.81. As shown in Paper I, the dependence of V^oo on 
(Si/s e )y/m e /mi is relatively weak. For this reason, we simply 
assume = 2.81 in this paper. 

In a weakly ionized gas, ^ generally depends on the ion- 
ization rate ( and the size distribution of dust particles. The 
value of *f> under a give condition can be cal culated from an 
algebraic equation (O09). We have shown in Paper I that the 
solution to the equation can be well fit by 



1 + 



-1/0.8 



(6) 



Here, 9 is a dimensionless quantity given by (IO09I : IPaper J) 

(n g e 2 



0: 
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where n g = p g / m g is the number density of gas molecules, and 
A lot and C tot is the total projected area and total radius of dust 
aggregates within a unit volume, respectively (note that we 
have assumed s,- = 1). In general, decreases with decreasing 
( and increasing the amount of dust. As found from Equa- 
tion ©, ^ approaches to only in the limit ^> and 
decreases with in the opposite limit. This reflects the fact 
that gas ionization is insufficient fo r dus t to be fully charged 
when is small (1O09I) . Following O09, we will refer to the 
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gas-dust mixture as the ion-electron plasma (IEP) and the 
ion-dust plasma (IDP) when > and O < ^>oo, respec- 
tively. 

If dust aggregates are monodisperse, we can write A tot = 
Auq/N and C tot = ano/N, where n ( > is the number density of 
dust monomers (note that uq/N represents the number den- 
sity of aggregates). Using these expressions as well as the 
scaling A = A/wa 2 ) and a = a^N 1 ! = aoN 1 / 2 described above, 
Equation (0 reduces to 



~MN) 



where 



Cn g e 2 
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(8) 



(9) 



is a dimensionless quantity representing the ionization rate of 
the gas. 

The electrostatic energy for monodi sperse a ggregates is 
given by (see Equations (27) and (29) of Pa per J) 



E E 



2 



N 



1/2 



where 



h = 



_ ^aoksT 



(10) 
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is the electrostatic energy between monomers in the IEP state 
(normalized by k%T). If we use Equations © and (8), Equa- 
tion ( [Tol l can be rewritten as 
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Thus, the normalized energies £o and £e are characterized 
by three dimensionless parameters (fo, h). 

2.2. The Drift Mass and Plasma Transition Mass 

It is useful to introduce two critical masses characterizing 
the motion and charge state of aggregates. The first one is 
the drift mass Md(= Nome,) defined by £o(No) = 1 (or equiv- 
alently, Eo{Mo) = k%T). This represents the mass at which 
the differential drift begins to take over Brownian motion in 
the collision velocity. Using Equation (fj), the equation for No 
can be rewritten as 



A(N D ) 2 



= foe 2 



(13) 



Thus, No is a function of foe 2 . As shown in PaperU No in- 
creases with decreasing foe 2 and behaves as No ~ l/b 2 foe 2 
in the limit of f D e 2 < 1, where b = 1/0.352 2.84 is the 
asymptotic value of N/A in the limit of N ^> 1 (see Equa- 
tion (HJi and the remark below the equation). Note that No ap- 
proaches unity as foe 2 — > 1 and becomes ill-defined at higher 
/be 2 . For this reason, we simply set No = 1 when foe 2 > 1. 

The second critical mass is the plasma transition mass 
Mp(= Npmo) defined by Q(Np) = ^oo- Since <d increases with 
N, the gas-dust mixture is in the IDP state when N -C iVp and 
is in the IEP state when N^> Np. With Equation ((H), the equa- 
tion for Np can be rewritten as 



A(N P ) 



N 



3/2 



= h. 



(14) 



TABLE 1 

Three Outcomes of the Growth of 
Charged Dust 



Conditions 




Outcome 


£e(N d )>6 




Total freezeout 


£e(N d )<6 **<*oo 


/4 


Bimodal growth 


£e(N d )<6 *»>*oo 


/4 


Unimodal growth 



The plasma transition mass increases with decreasing h and 
behaves as Np « 1 jb 2 h 2 in the limit of h <C 1 • Hence, if /V 
and Np ^> 1, we can approximate Equation ( fl2b as 
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We will use this approximation in the following subsection. 
Again, we set Np = 1 when h > 1 . 

2.3. The Growth Criteria 

Now we are ready to write down the gr owth cr iteria. The 
first criterion is given by (Equation (57) of iPaper t) 



£e(N d ) > 6. 



(16) 



If this inequality holds, all dust aggregates stop growing at a 
certain size ("total freezeout"). The size of the "frozen" ag- 
gregates is characterized by the freezeout mass Mp(= Nptno) 
defined by £e(Nf) = £k{Nf), where £k= 1 +£b is the total ki- 
netic en ergy (i.e., thermal energy + drift energy). As shown in 
IPaper 1 Mp is smaller than Mo whenever Equation ( TToT ) holds 
(i.e., the total freezeout occurs only in the Brownian motion 
regime), so the definition of Mp is effectively equivalent to 
£e(Nf)~L 

When the inequality in Equation ( TToT ) is reversed, the out- 
come of du st growt h depends on the second criterion (Equa- 
tion (58) of PapirJ 
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(17) 



If this inequality holds, a certain number of aggregates con- 
tinue growing while the rest stop growing at M sa Mo ("bi- 
modal growth"). Otherwise, all aggregates continue grow- 
ing with a relatively narrow size distribution ("unimodal 
growth"). The growth criteria are summarized in TableQ] 

For later convenience, we derive an approximate expres- 
sion of £e(Nd) applicable for foe 2 ,h <C 1. Recall that Nd(~ 
l/b 2 fo€- 2 ) and Np(~ l/b 2 h 2 ) are much larger than unity if 
foe 2 -C 1 and /i« 1, respectively (see Section 2.2). There- 
fore, we can approximate £ E (No) by Equation ( fT3T >. Thus, we 
obtain 



f 4N l J 2 . 



Se 



£e(N d ) : 
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No^>N P , 



Nd <C Np, 



(18) 



2be^f D 

for foe 2 , h <C 1. Equation ( fT8l is useful because £ E (No) is 
explicitly given as a function of (foe 2 , fE,h). 
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3. DISK MODEL 

The criteria shown in the previous section enables us to ex- 
amine how small dust aggregates evolve at each location in 
a protoplanetary disk. Here, we introduce a disk model used 
in this paper. The model is essentially the same as the one 
adopted in IO09L 

3.1. Structure of the Gas Disk 
We assume the gas surface density E g obeying a power law 



£,(/•)= 1.7 x 10V 



1 AU 



-3/2 



g cm 



(19) 



where r is the distance from the central star, and r/s is a scal- 
ing factor. The model with r^s = 1 is kn own as the mini mum- 
mass solar nebula (MMSN) model of lHavashil (11981b . We 
leave r/s as a free parameter to clarify the dependence on the 
disk mass. 

The gas temperature T is assumed to be isothermal in the 
ve rtical dir e ction, and the radial profile of T is given by that 
of lHavashil ( fl98ll) . 



7 = 280 



1 AU 



-1/2 _ 



K. 



(20) 



In reality, the temperature can be lower because of the large 
optical thickness provided by dust and can be higher because 
of the turbulent heating. However, we ignore these effects for 
simplicity. 

To obtain the vertical structure of the disk, we assume hy- 
drostatic equilibrium of the gas in the vertical direction. This 
leads to 



Pg : 



2nH 



exp 



2H 2 )' 



where z is the height from the disk midplane and 



H 



(21) 



(22) 



is the gas scale height. The isothermal sound velocity c s and 
the Keplerian orbital frequency fl^ are given by 




and 



GM* 



(23) 



(24) 



where G is the gravitational constant and M* is the mass of 
the central star. We assume a mean molecular weight of 2.34 
and write m g = 2.34toh, where toh is the hydrogen mass. For 
the stellar mass, we assume M* = 1M Q . 

Dust material is assume to be well mixed in the disk, and 
the dust density pd is related to p g as 



Pd(r,z) = fd g Pg(r,z), 



(25) 



where fy g is the dust-to-gas mass ratio. We choose fy g = 0.014 
as calculated from th e solar system abun dance of condensates 
including water ice (Poll ack et al.| [T994). The bulk density po 
of dust monomers is set to 1 .4 g cm -3 consistently with the 
adopted solar system abundance. We ignore the sublimation 
of water ice in inner disk regions for simplicity. 



3.2. Dust Motion 

We consider vertical sedimentation and disk turbulence as 
the mechanism driving differential drift between aggregates. 
Assuming z <C r, the vertical component of the stellar gravity 
is given by 



8s = M K z. 



(26) 



The relative velocity driven by turbulence generally depends 
on the ratio of the stopping times (r/,i and 7/2) of col- 
liding aggregates to th e turnover times of turbulent eddies 
(Orme l & Cuzzil 120071) . In particular, if the stopping times 
are shorter than the turnover time t v of the smallest eddies, the 
turbulence-driven relative speed can be approximately written 
in the form g T \ Tf t \ - tj^ | , where 



gT ~ Ur)/t v 



(27) 



is the effective acceleration driven by turbule nce, and u n is the 
characteristic velocity of the smallest eddies ( Weidenschilling 
IT984l:lOrmel & Cuzzill2007h . In this case, the total drift accel- 
eration g is given by 



2 2,2 
8 =8s + 8ti 



(28) 



where gs and gj are the vertical gravity and the effective ac- 
celeration driven by turbulence, respectively. The validity of 
the strong-coupling approximation (i.e., r/,i, 1/.2 <C t n ) is dis- 
cussed later. 

We evaluate gj in the following way. We express the 
strength of turbulence with the turbulent viscosity v tm \, = 
«Cj/57k, where a is the so-called alpha parameter for tur- 
bulence. The turbulent viscosity can be alternatively writ- 
ten as i4 ul b = u\ti, where ui and ti represent the character- 
istic velocity and correlation time of the largest eddies, re- 
spectively. We as sume t l = 1 /Qk as is for mag n etorotational 
turbu lence (e.g., [Fromang & Papa loizoul 120061; iTurner etldl 
2006). Equating the two expressions for z/ tur b, the character- 
istic velocity is obtained as ui = y/ac s . Assuming the Kol- 
mogorov spectrum, u 7] and t v are given in terms of u L and t L 



as u v = Re ' Ul and t v = Re ' ti, where Re = Vx m b/v mo \ is 
the Reynolds number. The molecular viscosity v mo \ is given 
by t'moi = u g /(2n s a mo \), where u g = ^%/ttc s is the molecular 
thermal speed and a m0 \ = 2 x 10~ 15 cm 2 is the m olecular col- 
lision cross section (Cha pman & Co wling 1970). Using these 
relations, we can rewrite the turbulence-driven acceleration 
gr (Equation d27b ) as 



g T w V^Re 1/4 c. s ft K . 



(29) 



Note that gT oc a 3 / 4 because Re oc a. In the following section, 
we treat a as a free parameter. 

Now let us check the validity of the strong-coupling approx- 
imation for the turbulence-driven relative velocity. It is useful 
to rewrite the Reynolds number as 



Re: 



aT, g a moi c -z-/2H 2 
2m g 



(30) 



where we have used v tm \, = ac s H, v mo \ = \f2/irc s m g /(p g <j mo \), 
and Equation d2Tb . Substituting this expression and Equa- 
tion (O into t v = Re" 1 / 2 ?,, = Re~ I//2 /f2 K , we have 



(31) 
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Note that Re (and hence J1k?^) is independent of T and 57k- 
For the stopping time t/, we use the fact that the mass-to- 
area ratio of fractal (D ~ 2) aggregates approaches to a con- 
stant in the limit N ^> 1 (see Section 2.2). Substituting Equa- 
tions (|2TT>— (|23l> into Equation dD and using N/A < b m 2.84, 
we obtain 



y ~ 2 £„ 



Js Vo.l urn/ Vl 



3/2 



//2H 2 



(32) 



,0.1 /xm/ VI AU/ 

Comparing Equation d32b with Equation OTb . we find that the 
strong coupling approximation is valid (i.e., r/ <C f,,) as long 
as a < 1(T 2 , fl < 1 Mm, and r < 100 AU. 

3.3. Gas Ionization 

For ioniz ing sources, we co ns ider Galactic cos- 
mic rays dUmebavashi & Nakanol [19811) . stellar X- 
ravs (|Igea & G lassgold 1999), and radionuclides 
dUmebavashi & Nakanol (20091) . Thus, we decompose 
the ionization rate as £ (cr + Cxr + Cra where £cr, Cxr, and 
(ra denote the rate of ionization by cosmic rays, X-rays, and 
radionuclides, respectively. We do not consider thermal ion- 
ization because it is negligible at gas te mperatures <C 10 3 K, 
or at heliocentric distances ^> 0.1 AU (lUmebay ashil 1 9 8 31) . 

The cosmic-ray ionization ra te is given by a fitting formula 
(Umebayashi & Nakano 2009) 



Ccr0",z) : 







CR.O 



-Et(r,z)/£cR 



+ e 



-[E s (r)-S;(r,z)]/E C R 



1 + 



1 + 



ScR 



- S y (r)-S^(r,z ) N | 

ScR 



(33) 

where Ccr.o = 1 -0 x 10~ 17 s" 1 is the cosmic-ray ionization rate 
in the interstellar space, £cr ~ 96 g cm" 2 is the attenuation 
length of the ionization rate, and 

p,(r,zVz'=^erfc(— J (34) 



is the vertical gas column density above altitude z. 
For the radionuclide ionization, we assume (ra 



7 x 



10 19 s l , which corresp onds to the ionization rate by a short- 
lived radionuclide 26 Al dUme bavashi & Nakano 200 9]) with 
the abundance ratio of 26 A1/ 27 A1 = 5 x 10" 5 dLee et al.ll 19771) . 
We neglect other short-lived radionuclides and all long-lived 
ones since they give only minor contributions. We also ne- 
glect the decrease in 26 Al due to the radioactive decay; this 
can be done as long as we consider an early stage of dust evo- 



lution within a timescale of 10 6 



yr. 



The stellar X-r a y ioni zation rate has been calculated by 
llgea & Glassgoldl Jl999) using the Monte Carlo radiative 
transfer code inclu ding Compton scatter ing. A useful fitting 
formula is given by iTurner & Sanol d2008h . 



U^Cxr.o^) 2 ( 



AU 

^ e -£j(r,z)/£xR 



+ e 



2 x 10 30 erg s" 1 

[£ s (r)-£+( r , z )]/ExR 



(35) 




heliocentric distance r [AU] 

FIG . 2 . — Radial profile of the ionization rate at z = H for tjy: = 1 . The thick 
solid curve shows the total ionization rate £ , while the thin solid, dashed, and 
dotted curves show the contribution from Galactic cosmic rays (CcrX stellar 
X-rays (Cxr), and radionuclides (Cra)- 

where Lxr is the X-ray luminosity of the central star, and 
Cxr.o = 2.6 x 10" 15 s" 1 and S XR = 8.0 g cm" 2 are the fitting 
parameters. This fitting f ormula approxima t ely re produces 
the k B T XR = 5 keV result of llgea & Glassgoldl d 19991) at > 

1 g cm" 2 where scattered hard (> 5 keV) X-rays are respon- 
sible for the ionization. At higher altitudes, Equation d35l > 
underestimates the ionization rate since it ignores the contri- 
bution of softer X-rays. We nevertheless use Equation ( 135] ) in 
this paper because the critical energy £e(Nd) is independent 
of £ at such high altitudes (see Equation d44li). We take Lxr = 

2 x 10 30 erg s" 1 in accordance with the median characteris- 
tic X-ray luminosity observed by Cha ndra for young s olar- 
mass stars in the Orion Nebula Cluster dWolk et alj2005l) . Al- 
though the c haracteristic X-ray temperature &b7xr ~ 2.4 keV 
observed by IWolk et ail d2005l) is lower than the assumed 
value of 5 keV, the choice of the temperatu re does not signif- 
icantl y affect the resulting ionization rate dlgea & Glassgoldl 

tun. 

As an example, Figure [2] shows the radial profile of the 
ionization rate measured at z = H for 77s = 1 . Cosmic rays 
dominate ionization at r > 20 AU, while X-rays dominate at 
4 AU < r < 20 AU. At r < 1 AU, both cosmic rays and X- 
rays are significantly attenuated, and thus the ionization rate 
reaches the floor value Cra- 

3.4. Global Profiles of fo, /e, and h 

With the disk model described above, we can obtain ana- 
lytical expressions for fo, /e, and h as a function of r and 

z. 

First, we decompose fo (Equation (0) into two compo- 
nents, 

!d= fD.S + fD.T, (36) 

where fo.s and fo.r are the contributions from stellar gravity 
an turbulence (i.e., fo with g = gs and g = gj), respectively. 
Using Equations (F2ll i-(l23l and d26i i. we can rewrite fo t s as 

97r 2 mo 



_ mp / mp y(Z_y z 2 ih 2 
lD - S ~ mmAnal^) \HJ & 



<Ax 10 



-i( «o \ 5 ( 
V0.1 urn) V 



/im/ V 10 3 g cm 



JC37) 



where the radial profile of Spis given by Equation d!9l . For 
fo.r we use Equations (f2TT > — (|23T>. (1291 . and (f30b to obtain 



97r 2 m r—( m 
/d,t = ™ — aVRel — j— 

128 m„ \-Kaf-Xi. 



z 2 /H 2 
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heliocentric distance r [AU] 

FIG . 3. — Radial profiles offo,s (Equa tion {57); solid black line), /e (Equa- 
tion (39); dashed line), and h (Equation )40t ; dotted line) atz = H for rj^ = 1 
and ao = 0. 1 /im. The two gray lines show foj (Equation 1381 ) for a = 1CT 4 
and 1(T 2 . 



VlO-V lo.l/mJ Vl0 3 gcm-V 



-3/2 



3z 2 /4tf 2 



(38) 



Note that both fo s and fo j are independent of T and J7k for 
fixed z/H. 
For Equation (fTTT i directly gives 



14 



0.1 /W V300K 



(39) 



where the temperature profile is given by Eq uation ( 1201 ). Fi- 
nally, for h, substitution of Equations (I2i1i-(|25T> and «o = 
pd/ffiQ into Equation (0 leads to 




e 



c 



■e 



: /2H 2 



/im/ V300K 

27r/»K ^ e ^/2H 

1 yr 



) ( 



10 3 g cm 2 - 



(40) 



As an example, we plot in Figure [3] the radial profiles of 
/d.s, /d,t, /e, and h&tz = H for = 1 and ao = 0.1 /im. It is 
seen that the ratio of foj / /d,s decreases with r, meaning that 
the effect of turbulent is relatively insignificant at outer parts 
of the disk. This is because the ratio fo,r / '/d.s is proportional 
to Re 1 / 2 and the Reynolds number Re cx E g decreases with r 
for fixed a and z/H. 

4. RESULTS 

4. 1 . Fiducial Case 

As the fiducial example, we begin with the case of t/s = 1 
(i.e., the original MMSN), a = (i.e., laminar disk), and ao = 
0.1 /im. 

The upper panel of Figure H] shows the radial profiles of 
£e(Nd) and calculated at one scale height above the mid- 
plane (z = H). We find that £e(Nd) exceeds the critical value 6 
at 1 AU < r < 100 AU. The total freezeout occurs in this wide 
region. The outcome of dust growth outside the frozen zone 



9? 10' 
^ 10° 

5 10- 2 

10" 2 
10" 4 
10~ 6 





y ^£e(W d )/6 












heliocentric distance r [AU] 

FIG. 4. — Upper panel: radial profile of £e(Nd) (solid curve) and Vf* 
(dashed curve) at z, = H for the fiducial model. At heliocentric distances 
where Se(Nd) 6 (shaded area), fractal dust growth stalls at mass M m Mf 
because of the electrostatic barrier. Lower panel: freezeout mass Mf (thick 
solid line), drift mass Mo (thin solid line), plasma transition mass Mp (dashed 
line), and critical restructuring mass M ro n (dot-dashed line) at z = H for the 
same model. The initial (monomer) mass mq is shown by the dotted line. The 
gray arrows illustrate the growth history of dust growth starting from vari- 
ous locations. Aggregates grown beyond M = M m \\ leave the fractal growth 
regime and will finally drift inwards as schematically shown by the dashed 
curved arrow (see also Section 5.2.2). 



depends on the value of (Equation (fT7l)). As seen in the fig- 
ure, increases with r and exceeds the critical value ^oo/A 
at r « 50 AU. Hence, the growth is unimodal at r > 100 AU 
and is bimodal at r < 1AU. 

The lower panel of Figure [4] plots the freezeout mass Mf 
in the frozen zone as well as the drift mass Mo and the 
plasma transition mass Mp at z = H. It is seen that Mf 
rapidly decreases as r increases, with the maximum value 
Mf ~ 10~ 7 g (a ~ 0.3 mm) and the minimum value as low 
as a few monomer masses. 

To see how early dust evolution stalls in the frozen zone, 
we compare Mf with the threshold mass M ro u for the onset 
of collisional compaction. The threshold mass is defined by 
E K (M m \\) = Eton, where 



Emu = 3tt 7an£crit 
«6xl0-'°( 



7 



100 erg cm 



/ Ccrit \ / «o 

\2 A/ V0.1 /im 



erg (41) 



is the energy needed to roll one monomer on another in con- 
tact by 90°, 7 is the surface adhesion energy for the two 
monomers, and £ CI it i s the critical tangential displacem ent for 
starting the rolling (jDominik & Tielensl [l 9951 1 1 9971) . Nu- 
meric al simulations (Dominik & Tielen slll997USuvama et alj 
2008) and laboratory experiments dBlum & Wurml |2000T) 
show that dust grows into fractal aggregates as long as the 
collision energy E K is below E m \\. Therefore, we can regard 
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?te=1, a =0.1yum, a=0 




heliocentric distance r [AU] 

FIG. 5. — Two-dimensional map showing the location of the frozen zone in 
the fiducial disk mode l. Shaded is the "frozen" zone, where the total freezeout 
occurs (Equation H6\ ). In this region, fractal dust growth stalls at the freeze- 
out mass Mf (whos e va lue is shown by the contours). The second growth 
criterion (Equation is satisfied above the dashed curve, indicating that 
the unimodal growth occurs outside above the frozen zone while the bimodal 
growth occurs inside the frozen zone. Outside the dotted curve, the drift mass 
Mo is larger than the p lasm a transition mass Mp, so the approximate formula 
for £e(Nd) (Equation d42» is applicable. 

M ro n as the maximum mass below which the collisional com- 
paction can be n eglected. For icy m onomers, 7 is estimated 
as 100 erg cm" 2 dlsraelachvilil 11992b but a realistic value of 
£crit is unknown. For a conservative estimation, we assume 
the minimum dis placement £ cr j t = 2 A an ticipated by the mi- 
croscopic theory (Domini k & Tielensl 19971) . which makes our 
aggregates the most easily compressed. In the lower panel of 
Figure @] we plot M ro ii as a function of r. We see that Mf is 
at least four orders of magnitude smaller than M ro u- This is 
a robust result because Brownian motion dominates the rela- 
tive motion of frozen aggregates and because the thermal en- 
ergy ~ kftT is generally much smaller than Ziroii- Note that the 
electrostatic barrier is in marked contrast to the other known 
growth barriers (e.g., radial drift, fragmentation, bouncing) 
which obstruct dust growth after the collisional compaction 
becomes effective. 

It is worth mentioning that the frozen zone always has in- 
ner and outer boundaries at finite r. In outer regions where 
Md 3> Mp (see FigureE}, the approximate formula for £e(Nd) 
(Equation ( fT8l ) reads 




(42) 

where we have used Equations (f37T > and (|39i l. Note that 
£e(Nd) at fixed z/H decreases as r increases because S ? and 
T generally decrease with increasing r. On the other hand, in 
inner regions where M D <C M P , £ E (No) is smaller than that 
in Equation g2) by a factor M D /M P oc CV^i^ 2 oc ( 2 r 3 /T 2 
(see Equation (fT8l). In this region, £ e {Nd) oc £ ? £ 2 r 3 /r de- 
creases as r decreases (unless S g /T is steeper than r" 3 ), since 
the ionization rate £ generally decreases as the surface den- 
sity increases. From these facts, we find that the frozen zone 
always has inner and outer boundaries at finite r. 

The two-dimensional (r— z) map of the frozen zone is dis- 
played in Figure |5] At z > H, the frozen zone moves inwards 



/fe=1, a =0.1jum, a=0 




heliocentric distance r [AU] 

FIG. 6. — Gas column density S glow of the dust growth zones at different 
heliocentric distances r for the fiducial model. The dotted line show the total 
gas column density S g . 

as z increases, because £e{Nd) is lower at larger z and r (see 
Equation d42i>). A region very close to the midplane is en- 
tirely covered by the frozen zone because the settling veloc- 
ity vanishes there (i.e., £jy — > 0). The bimodal growth zone 
covers middle altitudes at r < 1 AU and expands towards 
smaller r. The upper boundary between the bimodal and 
frozen zones is roughly characterized by the column depth 
£g(z) ~ 10 2 g cm" 2 . This reflects the fact that cosmic rays is 
significantly attenuated at such depths. The unimodal growth 
is only allowed at large r and z where the settling velocity is 
high enough for dust to overcome the electrostatic barrier. 

It will be useful to show what amount of dust is allowed 
to grow at different heliocentric distances. Figure [6] plots the 
column density £ grow of gas within the growth zones (bimodal 
and unimodal zones) as a function of r for the fiducial model. 
Note that the column density of dust within the growth zones 
is /dgSgrow, while the column density of gas within the frozen 
zone is £ g — E grow . For comparison, the total gas surface den- 
sity S ? is also shown by the dotted line. We see that £ grow 
is comparable to £„ at r < 1 AU because of the presence 
of a large bimodal zone (see Figure [5}. This means that the 
electrostatic barrier does not strongly affect dust growth at 
r < 1 AU. Farther out the disk, however, E grow steeply de- 
clines to lO" 4 !^ fa 10" 1 g cm . The residual value of E„ rnw 
comes from the unimodal zone at high altitudes. At 1-10 AU, 
the amount of dust within the growth zone is less than 1% of 
the total dust mass. 

Figure [6] implies that the column density of the unimodal 
zone depends on r only weakly. As shown below, this comes 
from the weak radial dependence of the gas temperature T. 
First, we note that the gas column density above altitude z, 
S*(z) (Equation d34li), approximately behaves as 

£; ( z)«^e-/ 2 " 2 (43) 
s V2tt z 

at z 3> H (this equation directly follows from the asymptotic 
expansion of the complementary error function erfc(x)). Us- 
ing this approximation, Equation (|42| > can be rewritten as 

^)«35(-^)- 3/2 (-l-)(-^). (44) 
V0.1 jitm/ V 100 K/ V 1 g cm 1 ) 

This equation means that, for fixed r, the column density 
above a given altitude z is proportional to £e(Nd) at that al- 
titude. Substituting £e(Nd) = 6 and = E grow /2 into this 
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FIG. 7. — Same as Figure[6] but for different values of disk mass. The black 
solid and dotted curves indicate S grow and T, g for rfs = 0. 1 , respectively. The 
gray curves are for the fiducial model of rjs = 1 ■ 
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FIG. 8. — Same as Figure [6] but for different monomer radii ag. The 
black solid and dashed curves show Egrow for uq = 0.5 fim and 0.02 fim, 
respectively. The gray curve is for the fiducial model of ao = 0. 1 fim. 
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equation, we obtain the column density E er( 
growth region as 

3/2 , J 

V100K 



of the unimodal 



-"grow 1 



V0.1 //m 



gem 



(45) 



Thus, we find that the column density of the unimodal zone 
is inversely proportional to T and hence depends on r only 
weakly. It should be noted that the column density of the 
unimodal zone is independent of Q, 'Eg, and fd g . 

4.2. Dependence on the Model Parameters 
4.2.1. Disk Mass 

To examine the effect of the disk mass, in Figure|7] we plot 
the radial profile of S glow for a lighter disk model of r)s = 0. 1 . 
As in the fiducial case, the high value of £ grow at small r corre- 
sponds to the bimodal zone, while the floor of S grow at larger 
r comes from the unimodal zone. We find that the bimodal 
zone shrinks as the disk mass decreases. This is because ion- 
ization sources from upper layers (cosmic rays and X-rays) 
can penetrate to the midplane more easily when the disk mass 
is smaller. However, this effect is limited because a smaller 
disk mass also causes a higher settling velocity and hence a 
higher collisional velocity. Moreover, the column density of 
the unimodal zone is very insensitive to 77s as is already ex- 
pected from Equation d45b . Thus, we can conclude that the 
change in the disk mass has only a minor effect on the size of 
the frozen zone. 

4.2.2. Monomer Size 

We have assumed that dust growth begin with single-sized 
monomers. In reality, however, dust monomers in protoplan- 
etary disks will obey a certain size distribution. For instance, 
interstellar extinction i mplies that the gra in sizes range from 
0.005 /zm to 0.25 ztm (iMafhis et al.lll977l) . To fully take into 
account the monomer size distribution is challenging since 
the porosity model for aggregates would be much more com- 
plicated. For this reason, we try to estimate the effect of 
monomer size distribution by varying the monomer size ao 
in our model. 

Figure [8] shows £ grow for three different monomer radii, 
ao = 0.02, 0.1, and 0.5 /im. We see that the increase of ao 
by a factor of 5 leads to the enhancement of the column den- 
sity of the unimodal zone (i.e., E grow at r > 1 AU) by a fac- 
tor of 10. This is consistent with Equation (|43T > predicting 
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FIG. 9. — Same as Figure|5] but for different turbulence strengths a. The 
upper and lower panels are for a = 10~ 3 and 10~ 2 , respectively. 



3/2 

that the column density scales as a . The growth zone in- 
creases with fl because aggregates made of larger monomers 
have larger mass-to-area ratios and hence higher drift veloc- 
ities. However, even with ao = 0.5 /im (which is twice the 
maximum size in the MRN distribution), the enhanced S grow 
is still much smaller than the total surface density S g . Indeed, 
ao must be as large as 10 /im for S grow to be comparable to 
E g at every r. Thus, a change of the monomer size within a 
realistic range gives only a minor effect. 

4.2.3. Turbulence 
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Same as Figure|6] butfor different turbulence strengths a. The 
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black solid and dashed curves correspond to a = 10~ 2 and 10 
The gray curve is for even weaker turbulence, a ^ 10 



, respectively. 



Turbulence enhances the collision velocity of aggregates, 
so it acts to reduce the size of the frozen zone. To see this 
effect, we have computed S glow with changing the value of a. 
The two panels in Figure|9]sriow the two-dimensional maps of 
the frozen zone for a = 10~ 3 and 10~ 2 , and the radial profiles 
of E grow for the two cases are plotted in Figure [10] We have 
also computed Sgrow for a ^ 10~ 4 , but the result is indistin- 
guishable from that for the fiducial laminar case. Comparing 
Figure [9] with Figure [5] we find that the presence of turbu- 
lence is particularly important near the midplane. This is be- 
cause the settling velocity vanishes at the midplane while the 
turbulence-driven velocity does not. For a = 10~ 3 and 10~ 2 , 
the bimodal zone expands out to 3 AU and 10 AU, respec- 
tively. This means that strong turbulence with a ~ 10~ 2 re- 
moves the frozen zone from planet-forming regions. At larger 
r, however, removal of the frozen zone requires even higher 
values of a. We find that a must be higher than 10 _1 for E grow 
to be comparable to £„ at every r. 

Unfortunately, it is uncertain what is a realistic value of a 
for protoplanetary disks, especially in early stages of dust evo- 
lution. It is commonly believed that the most promising mech- 
anism drivi ng disk turbulence is the magnetorotational insta- 
bility (MRF lBalbus & Hawlevlll991l) . MHD simulations sug- 
gest that the MRI can sustain turbulence at a level of a ~ 10~ 2 
or hig her depending on t he net vertical flux of the magnetic 
field (ISuzukietaLlbOlOl) . This implies that the MRI could 
assist dust to overcome the electrostatic barrier at r < 10 AU 
(see Figure [Toll. The problem is whether the MRI does oper- 
ate there in early stages of dust evolution. Since the MRI is 
an MHD phenomenon, it requires an ionization degree high 
enough for the gas to couple to the magnetic fields. The MRI 
is suppressed by Ohmic dissipation where the ionization de- 
gree is too low, which is known as the "dead" zone. Im- 
portantly, the size of the dead zone can be very large when 
small and fluffy dust aggregates are abundant because they 
have large surface areas and therefore efficiently capture the 
ionized gas. For example, if all dust parti cles are 0.1 /jm i n 
size, the dead zone can extend to 20 AU dSano et al.l l2000). 
Furthermore, the large dead zone remains while the particles 
grow into fractal aggregates because their total projected area 
is nearly conserved (O09). These facts suggest that MRI- 
driven turbulence may be absent from inside 10 AU in the 
early stage of dust evolution. However, the physics deter- 
mining the size of the dead zone is not yet fully understood 
(see, e.g., Inutsuka & Sano 2005), so further investigation is 
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FIG. 1 1 . — Gas column density £ gm w of the growth zones as a function 
of r for the compact growth model (solid curve). For comparison, the result 
for the fiducial porous model is indicated by the grey curve. The dotted line 
show the total gas column density E s . 



needed on this issue. 

5. DISCUSSION 
5.1. On the Role of Porosity Evolution 

So far, we have assumed that dust grows into fractal ag- 
gregates of D ~ 2. The assumption of fractal evolution is 
reasonable when analyzing the electrostatic barrier since the 
freezeout occurs much earlier than the onset of collisional 
compaction (see Section 4.1). However, one cannot rule out 
that some unknown process could cause compaction of aggre- 
gates. Therefore, it will be useful to see how the result in the 
previous section depends on the porosity model adopted. 

Here, we consider the conventional, compact growth 
model where aggregates grow as zero-porosity spheres (e.g . , 
iNakagawa et alll98UlTanaka et al.l2005HBrauer et al. 2008). 
The effect of the electrostatic barrier on the compact dust 
growth has been already investigated in our previous paper 
(Paper I). We have found that the freezeout criterion for the 
compact model can be again written as Equation ( [T6l l, but, 
the electros tatic en ergy £e{Nd) is now given by (see Equa- 
tion (66) of iPaper t) 



£ E (N D ) = 



Se 



i+(hf D 3/ Y s 



-2.5 



1/5 



(46) 



where the definitions of fo, /e, and h are the same as those 
for the porous case (i.e., Equations (0, dTTT >. and ©, respec- 
tively). It can be easily checked from Equations (f37T> — (f40l> 
that the right-hand side of Equation d46l > is independent of the 
monomer radius oq. This should be so since "monomers" are 
not well defined in the compact model. 

Using the freezeout criterion and Equation (l46b . we can 
compute the column density £ grow of the growth zone in the 
same way as we did in Section 4. Figure [TTJ plots E grow for 
the compact dust model with the fiducial parameter rjs = 1 and 
a = 0. For comparison, we overplot the result for the fiducial 
porous model shown in Figure [6] We see that the frozen zone 
predicted by the compact model is considerable smaller than 
that by the porous model. This reflects the fact that the dif- 
ferential drift velocity of compact "aggregates" grows more 
rapidly with N and h ence takes over the thermal velocity at 
lower N (see iPaperll) . 

From the above example, we can say that the size of the 
frozen zone is quite sensitive to porosity evolution. However, 
it should be kept in mind that evolution into highly porous 
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aggregates are ine vitable unless some compaction mecha- 
nism exists (e.g. [Wurm & Bluml 119981: iBlum et all 1 19981: 
iKempf et al.lll999t IOTS09I) and that collisional compaction 
is not the mechanism in the early stage of dust evolution. Of 
course, there remains a possibility that there exist any other 
effective mechanisms. For example, compaction of fluffy ag- 
gregates is in principle possible if a number of small particles 
are supplied by some process and fill in the voids in the ag- 
gregates. To specify such mechanisms is beyond the scope of 
this paper, so we leave this issue open for future work. 

5.2. Effects of Global Dust Transport 

In Section 4.2.3, we have seen that turbulence as strong as 
a > 10~ 2 is preferable for dust growth beyond the electrostatic 
barrier. However, from the viewpoint of planetesimal for- 
mation, strong turbulence is not always preferable because it 
causes the fragmentation barrier against large and compacted 
aggregates. For example, large and compacted aggregates 
with Tf ~ fj" 1 acquires a random velocity of Am w m l ss ^/ac s 
in turbulence (e.g. JOrmel & Cuzzill2007l) . If we put a > 10" 2 
and c s « 700 m s" 1 (as is for T s» 100 K), the random veloc- 
ity amounts to Am > 70 m s -1 . By contrast, it is suggested 
from computer and laboratory collision experiments that the 
threshold collision velocity for catastrophic disruption to oc- 
cur is « 50 m s" 1 for icy aggregates ( Wada et al. 20091) and 
is an order of magnitude even lower for rocky aggregates 
dBlum & Wurm 112008b iGtittler et al.ll2010l) . Thus, even if we 
assume a high threshold velocity of ss 50 m s" 1 , the large ag- 
gregates cannot grow beyond the fragmentation barrier. 

For this reason, it is worth considering whether planetes- 
imal formation is possible in the frozen zone without strong 
turbulence. It is clearly difficult to form planetesimals directly 
from fluffy frozen aggregates, because fluffy aggregates cou- 
ple to the gas too strongly to form gravitationally unstable 
dust layers or clumps. Here, we examine a possibility that 
large and compacted aggregates are supplied from outside the 
frozen zone through some dust transport mechanism. 

For such transport mechanisms, we focus on (1) vertical 
mixing of frozen aggregates by (weak) turbulence and (2) ra- 
dial infall of aggregates grown at outer (> 100 AU) regions. 
As shown below, both mechanisms can supply large and com- 
pact aggregates to the frozen zone on a timescale of 10 6 yr or 
longer. 

5.2.1. Vertical Mixing by Weak Turbulence 

As seen in the previous sections, turbulence with a < 10~ 4 
does not significantly contribute to the collision velocity of 
aggregates. However, turbulence does not contribute only to 
the collisional velocity but also to the vertical mixing of dust 
material. Vertical mixing is efficient even by weak turbulence 
because fluffy aggregates strongly couple to the gas. The mix- 
ing will allow the frozen aggregates to go out of the frozen 
zone and grow there until reentering the frozen zone (see Fig- 
ure [T2l . This cycle can prevent the aggregates from being 
perfectly frozen. 

To fully examine this effect, we would have to solve the co- 
agulation equat ion including vertic al diffusion (as done by, 
e.g., |f)ullemond & Dominik 2005), but this is beyond the 
scope of the present study. In this paper, we restrict our- 
selves to simply estimating how long time is required for dust 
in the frozen zone to grow beyond the fractal growth regime 
(M > M ro \\) with the help of the turbulent diffusion. To do so, 
we estimate the mean growth time (w mean collision time) of 



turbulent mixing 



(^star^) 



disk 

growth & 
o °o compaction 

V settling 



'frozen" zone 



.O 



O 



-1 AU 



radial infall 



-100 AU 

FIG. 12. — Dust transport mechanisms that could supply large aggre- 
gates to the frozen zone (not to scale). Without any transport mechanism, 
the electrostatic barrier halts dust growth at the fractal growth stage at 1— 
10 AU < r < 100 AU near the midplane (the "frozen" zone). However, verti- 
cal mixing due to turbulence could allow the frozen aggregates to grow out- 
side the frozen zone. Furthermore, already large and compacted aggregates 
can be supplied from outer (> 100 AU) regions to the frozen zone as a re- 
sult of their radial infall. Both processes can provide large and compacted 
aggregates on a timescale of > 10 s yr (see text). 



aggregates at fixed r, 



' grow 



1 dM\ 
~M~dt) 



(47) 



where M is the mass of the aggregate. If the electrostatic 
barrier were absent, we would have dM/dt = pdcr co \\/S.u, and 
hence r grow = M / pda co nAu, where pd is the dust density, Am 
is the collision velocity, and cr co u is the collisional cross sec- 
tion. Here, we need to take into account the fact that aggregate 
collision is only allowed in the growth zone with the gas col- 
umn density £ grow . Assuming that the dust is vertically well 
mixed by turbulence, the probability that an aggregate experi- 
ences a collision in the growth zone will be given by the ratio 
Sg/Sgrow Therefore, we evaluate r grow as 

M 



LIH'U . . ' * (48) 

PrfCr co iiI\U Z^grow 

To simplify Equation d48l ). we use the fact that aggregates 
grow mainly through similar-sized collisions (OTSOy) and 
approximate ct co u as 47ra 2 . A collision takes the longest 
time when the relative velocity Am minimizes. Neglecting 
the contribution of weak turbulence to Am, the minimum 
value of Am is evaluated as the differential settling velocity, 
Am s» egTf s» eO|Hr/. Using r/ « M/Y> g n K A, A s» ira 2 , and 
p d ~ fdg^ g /H, we obtain 



1 



avow 

•10- 1 



;4x 10 H 



3/2 £„ 



100 AU, 



yr. (49) 



Note that the right-hand side of Equation (l49b does not involve 
M. 

In Figure [13] we plot r„ mv , as a function of r for a weakly 
turbulent case of a $J 10~ 4 . The jump in r erow at r s» 1 AU 
corresponds to the outer edge of the bimodal growth zone (see 
Figure [Toll. At 1 AU < r < 100 AU, where E grow < £ g , the 
typical value of r grow is about 10 5 yr (at 100 AU) to 10 6 yr (at 
1AU). 
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FIG. 13. — Timescale of dust growth driven by turbulent diffusion, Tg r0 w 
(Equation {49); thick solid curve). The turbulence parameter is set to a ^ 
10" 4 . The thin solid curve shows t„ iow when the frozen zone is absent. The 
parameters rfs and oq are the same as those for the fiducial model. 



Note that r grow only represents the mean time spent for a 
single collision. What we really want to know is the time re- 
quired to grow beyond the fractal growth regime. Using the 
fact that r grow = dt/d\\\M is independent of M at M > M D 
and is shorter at M < Mo, the required time is estimated as 
Tgrow hi(Mroii /Mo)- As found from the lower panel of Figure|H 
this time is longer than r grow by a factor of ~ 10. Hence, we 
conclude that turbulent mixing could assist the frozen aggre- 
gates to grow beyond the fractal growth regime on a timescale 
of>10 6 yr. 

It is interesting to note that the vertical mixing predicts 
a quite fiat radial profile for the dust growth timescale at 
1 AU < r < 100 AU. In fact, it can be directly shown from 
Equations d45b and d49l that the growth time at the locations 
is proportional to E ? /f2KSg row oc TY, g /Q^ and hence scales 

with r -1 / 2 in the MMSN modefl The flat profile is in marked 
contrast to previous uncharged dust models where r grow scales 

with the orbital period (Sl^ 1 ) and is hence as steep as r 3 ' 2 . 



by 



Tinfall = 




Dust aggregates continue growing at the same location until 
Tinfaii becomes shorter than the local growth timescale T grow 
(see Section 5.2.1 for its definition). Hence, the lower limit 
of u r for the drifting aggregates can be evaluated from the 
inequality T in f ail < r grow . At r w 100 AU, this inequality gives 
u r >5m s -1 . Thus, the lower limit of £k is estimated as 

Ek > \m f u] > lO" 9 (j^) erg, (51) 

where we have used the fact that the mass Mf of the frozen ag- 
gregates is much smaller than that of the drifting aggregates. 
For the electrostatic energy, we can obtain its upper limit by 
setting Q ^ ^ ooak^T / e. Assuming that the frozen aggregates 
is much smaller than the drifting aggregates, we obtain 

E E <a F [ J 

~'°-"(iif(4) 2 - <*> 

where a F = ao(Mp /mo) 1 / 2 is the radius of the frozen aggre- 
gates. Since Ek is much higher than E F , we conclude that 
the radially drifting aggregates can overcome the electrostatic 
barrier to collide with the frozen aggregates. 

However, it is unclear from the above argument whether the 
collision results in the growth or fragmentation (erosion) of 
the drift ing aggregates. Interestingly, recent labora tory exper- 
iments (iTeiser & Wurml 120091: IGuttler et al. 2010) show that 
net growth of a large and compact aggregate is possible even 
at high collision velocities (> 1 m s" 1 for rocky aggregates) if 
the projectiles are muc h smaller than the t arget (see the "pC"- 
regime in Figure 11 of IGuttler et a l. 2010). Hence, the pres- 
ence of frozen aggregates could even assist the growth of the 
drifting aggregates beyond the fragmentation barrier. 



5.2.2. Radial Infall of Large Aggregates from Outer Regions 

Aggregates outside the frozen zone are allowed to grow be- 
yond M = M ro \\ and experience collisional compaction. As 
the compaction proceeds, they will gradually decouple from 
the g as, settle onto the mid plane, and drift towards the central 
star (fWeidenschilhn gl fl 9771) . Aggregates falling from outer 
(r> 100 AU) growth zones should enter the frozen zone (see 
Figure [T2l. The time spent for the onset of the radial drift will 
be comparable to that for the onset of the collisional com- 
paction, which is ~ 10r glow ~ 10 6 yr at r w 100 AU (see Sec- 
tion 5.2.1). Thus, large and compact aggregates can be sup- 
plied from the outer regions on a timescale of ~ 10 6 yr. 

Furthermore, the drifting aggregates have a sufficiently high 
kinetic energy to sweep up frozen ones beyond the electro- 
static barrier. To show this, we compare the impact energy Ek 
with the electrostatic energy E F for the collision. The lower 
limit of Ek can be estimated by considering the lower limit of 
the radial drift velocity u r for the drifting aggregates. Here, it 
is useful to note that the timescale of the radial infall is given 

5 Tg ro w is indepen dent of r if S g oc r~ l as suggested by subrrdll ime- 
ter observations (e.g., Andrews & Williams 2007; Andrews et al. 2009) and 
constant-a accretion disk models (e.g., Hartmann et al. 1998). 



5.2.3. Implications for the Timescale of Protoplanetary Dust 

Growth 

As discussed above, dust evolution in the frozen zone is 
possible on a timescale of 10 6 yr or longer. By contrast, if 
the electrostatic barrier is absent, the collision timescale for 
fractal aggregates at r < 10 AU is as short as < 10 3 yr (see 
Figure [131 . Thus, the electrostatic barrier can dramatically 
alter the timescale of dust evolution in protoplanetary disks. 

There is possible evidence that evolution of small 
grains/aggregates occurs on a long timescale. Infrared ob- 
servations of classical T Tauri stars show mid-infrared exces s 
evolving on a timescale of 10 6 yr (e.g.. iFurlan et al.l 12006). 
This is usually interpreted as an indication that micron-sized 
and warm (> 100 K) gra ins grow in their curcumstell ar disks 
on this timescale (e.g., Dullemond & Dominik 2005). How- 
ever, it is also possible to interpret this timescale as the du- 
ration of fractal dust growth because the spectral signature 
of large and fractal (D ~ 2) ag gregates is sim ilar to that 
of small and compact particles dMin et all l2006t). Interest- 
ingly, the presence of the frozen zone at r > 1 AU together 
with the dust transport mechanisms considered above can 
explain the retention of fractal aggregates in warm regions 
(< 10 AU) on a timescale ~ 10 6 yr. Previously, the retention 
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of fast-growing small grains has been attributed to the col- 
lisional destruction of large aggregates into fragments (e.g. , 
Dullemond & Dominikl l200lt iDominik & Dullemondl 120081: 
Birnst iel et al.1 120091) . However, the electrostatic barrier can 
be an alternative solution to this problem. 

Another piece of possible evidence for slow dust evolu- 
tion may be obtained from primitive meteorites. Radioiso- 
tope dating of chondrites has revealed that the the formation 
of chondrules began 1-3 mill ion years after t he for mation 
of Ca-Al-rich inclusions (e.g. iKita et all l2000[ l200l . This 
implies that the growth of chondrule precursors (i.e., mm- 
sized dust aggregates) occurred on a timescale of 10 6 yr. The 
electrostatic barrie r or the "bouncing barrier" proposed by 
IZsom et ail d2010l) might have played a role in this slow ac- 
cretion process. 

6. SUMMARY 

In this paper, we have examined where in a protoplane- 
tary disk dust charging can halt local dust evolution. This 
is the first step towards the modeling of dust evolution in 
protoplanetary disks including dust charging together with 
other importan t mechanisms, such as collisional comp action 
(ISuvama et al.l 120081) . radial drift (IWeiden schilling Il977h 
Braue ret al.l l2008|Kbouncing (|Zsom et al.l 120101) . and frag- 
mentation dBrauer et aT]2 008; Birn stiel et al.1 2009). Our find- 
ings are summarized as follows. 



1. We find a "frozen zone" where dust growth stalls in 
the fractal growth stage. For weakly turbulent disks 
(a < 10~ 2 ), the frozen zone contains the major part of 
dust materials at a few AU to 100 AU from the cen- 
tral star. Dust growth beyond the fractal stage is only 
allowed in an inner region where ionizing cosmic rays 
and X-rays do not reach, and in an outer region where 
settling velocity is high enough to overcome the elec- 
trostatic barrier. The freezeout mass, the mass at which 
the growth begins to stall, strongly depends on the 
distance from the central star, typically ranging from 
10" 7 g (at - 1 AU) to as small as 10" 13 g (at > 10 AU). 

2. The size of the frozen zone does not significantly 
change with changing the disk mass and monomer size 
within a realistic range (Sections 4.1 and 4.2). By con- 
trast, turbulence as strong as a > 10~ 2 can help dust in 
the "planet-forming" region (< 10 AU) to overcome the 



electrostatic barrier (Section 4.3). Caution is needed 
when considering this result in the context of planetes- 
imal formation, since such strong turbulence can cause 
the fragmentation barrier after the fractal growth stage. 

3. For weakly turbulent disks, we have considered two 
dust transport mechanisms that could lead to dust evo- 
lution in the frozen zone (Sections 5.2). Turbulent mix- 
ing across the boundary of the frozen zone can pre- 
vent dust growth from being completely frozen. Large 
and compacted aggregates can be also supplied from 
large heliocentric distances (> 100 AU) through their 
radial infall. Both mechanisms can result in the sup- 
ply of large and compacted aggregates and the re- 
moval of small and fractal aggregates on a timescale 
of 10 6 yr or longer. This is in clear contrast to previ- 
ous theoretical understanding that without fragmenta- 
tion, small dust particles get depleted in disks on much 
shorter timescales (e.g.. iDullemond & Dominikl 12005b 
iDominik & Dullemondll2008l) . This might explain the 
"slow" (~ 10 6 yr) dust evolution su ggested by infrared 
observation of T Tauri stars (e.g., iFurlan et al.l 2006) 
and b y radioactive dating of chondrules (e.g. jKita et alj 
120051) . 

Finally, we remark that the above findings are obtained us- 
ing a fractal growth model for small dust aggregates. This is 
in marked contrast to most previous studies on dust coagula- 
tion where aggregates are simplified as compact spheres (e.g ., 
iNakagawa et alJl98UlTanaka et al.l2005llBrauer et al.1 2008). 
However, it is critical to properly take into account porosity 
evolution because compact spheres would be more resistive 
to the electrostatic barrier (Section 5.1). On the other hand, it 
is also true that the effect of the electrostatic repulsion could 
be reduced if some mechanism prevents aggregates from be- 
ing highly fluffy (D ~ 2). As shown in Section 4.1, collisional 
compaction is very unlikely to be the mechanism because the 
freezeout begins much earlier than it becomes effective. Nev- 
ertheless, we cannot rule out possible existence of any other 
compaction mechanisms. We leave this issue open for future 
work. 

We thank the anonymous referee for careful reading of the 
manuscript and useful comments that helped improve it. S.O. 
is supported by Grants-in-Aid for JSPS Fellows (22 • 7006) 
from MEXT of Japan. 
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